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Abstract 

We theoretically study cooperative effects in the steady-state trans- 
mission of photons through a medium of N radiators. Using methods 
from quantum transport, we find a cross-over in scaling from N to N 2 in 
the current and to even higher powers of N in the higher cumulants of the 
photon counting statistics as a function of the tunable source occupation. 
The effect should be observable for atoms confined within a nano-cell with 
a pumped optical cavity as photon source. 

1 Introduction 

Collective effects in the emission of light have attracted a lot of attention re- 
cently. Super-radiant scattering off Bose-Einstein condensates [1], collective 
self-organization of atoms [2, 3], or the observation of the phase transition [4] 
for ultra-strong coupling in a cavity [5] are some of the modern manifestations 
of the original Dicke effect [6, 7]. 

There has also been an increased interest in identifying similar collective 
interference phenomena in mesoscopic (electronic) transport [8] . Many of these 
activities are triggered by progress in the time-resolved detection and counting 
of individual electrons [9] , and the fabrication of quantum circuits as electronic 
test-beds for quantum optical effects [10]. 

With the progress obtained in single-photon detectors [11] we feel that it is 
well motivated to revert the focus by applying theoretical methods from Full 
Counting Statistics (FCS) [12] to steady-state transmission of photons through 
a medium of N radiators in the small-sample limit of the original super-radiance 
model [13]. We show that higher cumulants of the photon counting statistics 
are sensitive indicators of cooperative emission effects even for only moderate 
pumping of the radiators. 

Steady-state super-radiance of incoherently pumped atoms in a cavity has 
been analyzed recently by Meiser and co-workers in terms of time-dependent 
intensity correlation functions [14]. Our results provide an alternative scheme, 
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where we propose to confine atoms within a nano-cell and to use a pumped op- 
tical cavity as photon source. The radiation emitted by the atoms is analyzed in 
terms of cumulants of a stationary photo-detection statistics. Instead of directly 
evaluating the cumulant-generating function - this is analytically impossible for 
the large system sizes considered - we extract the long-term cumulant behavior 
from the Laplace transform of the moment-generating function. We also expect 
our findings to be relevant for other situations dealing with collective boson 
transport, such as thermal transport [15] or phonon lasers [16]. 



2 Model 



We consider a number of N a bosonic reservoirs a connected via a medium of 
N two-level systems with identical level splittings f2. The collective coupling is 
described by an extended multi-mode Dicke Hamiltonian (cf. Figure 1 a), 
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with collective pseudo-spin operators J x ' z = Y^^=i <T i' Z ■ Here, b\ a creates a 
boson mode k with frequency wt Q in reservoir a and h^a denote the coupling 
constants. We derive the counting statistics of the total number n of bosons ex- 
changed between the medium and one (a) of the reservoirs by formally introduc- 
ing [17] a counting operator d\ = ^2n=-oo \ n + -0 ( n \ ( an( i similarly for the re- 
verse tunneling processes) in the Hamiltonian via b* k& —> b\-®d & . In all what fol- 
lows, we consider the weak coupling regime and derive a Lindblad master equa- 
tion in Born, Markov, and secular (BMS) approximation describing the time 
evolution of the system (TV two-level medium). The system state p^ = (n\ p \n) 
is conditioned upon the number of particles measured in the detector n. The 
master equation assumes the form p^ = Cop^ + £+p' n-1 ' + C~p( n+1 \ where 
jumps into (from) the bath a are described by the super-operator £ + (£_) and 
the probability of having counted n particles after time t is given by P n (t) = 
Tr { p( n ) (t) } . These n- resolved master equations are Fourier-transformed by in- 
troducing a counting field \ p(x) = P emx > an( i similar for further 
detectors possibly placed in the other baths, leading to a generalization of the 
master equation for collective spontaneous emission and re-absorption [18], 
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where we have omitted the dependence of p on the counting fields for brevity 
and r a = 27r^ fc \hk a \ 2 8(£l — ujka) denotes the spontaneous photon emission rate 
of a single two-level system into reservoir a, for which we assume stationary 
occupations n a = (b\ a bk a ) of photons with frequency ft = <jfc. In the case of 



thermal photon reservoirs, one has n a = [< 
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For a single counting field x> cumulants (denoted with double brackets 
throughout) of the photo-detection probability distribution P n (t) of n photons 
after time t are obtained as ((n fe (i))) = (— id x ) k C(\, t)\ x= o, where C(x,t) = 
hiTr {p(x,t)} = InTr {e^^'po} is the cumulant generating function. 



3 Results 

3.1 Single thermal reservoir 

For a single (N a = 1) thermal bath at inverse temperature (3, the Liouvil- 
lian C(x = 0) fulfills detailed balance and leads to equilibration of the sys- 
tem temperature with the stationary state of Eq.(2) is simply given by 
p oc e ~P^/ 2J , l n the angular momentum eigenbasis (J z \ jm) = 2m \ jm) with 
j = N/2 and -N/2 <m< +N/2), C 

has tridiagonal structure, and we obtain the FCS for transient Dicke super- 
radiance [13] with all TV atoms initially in the excited state. At large times t, 
the FCS is given by the probability of n photon emission and N — n absorp- 
tion events, i.e., the distribution (for < n < N) approaches P n (t — > oo) oc 
[1 + ns]™ n B~ n wrtn stationary moments (n fe ) = ^lim^ (— id x ) k Tr {p(x, t)} | x =o 



= T 1 =H N - ( J V2 + N/2)] k ) _ , (3) 



which explicitly shows that the statistics of emitted photons and the statistics 
of excitations remaining in the thermalized system are perfectly correlated. In 
a vacuum bath (/3f2 — > oo), all N initial excitations simply tunnel out of the 
system. For very large bath occupations (/3f2 —> 0) however, the correspond- 
ing higher cumulants grow super-lincarly with the system size ((n)) = N/2, 
((n 2 )) = N(N + 2)/12, ((n 4 )) = -N{N + 2) [N(N + 2) + 2] /120 and so on 
(all odd higher cumulants vanish). This super-linear scaling clearly is a col- 
lective effect: If the N two- level atoms decayed independently, all cumulants 
would scale linearly in N (cumulants of independent processes are additive). 
Measuring these higher cumulants would require to detect all emitted photons 
with high accuracy. 

For measurements of transient FCS at finite times t < oo, however, it is 
also necessary to consider the finite bandwidth of the photo-detector given by 
its time resolution At. The probability of detecting n photons during the 
time interval [t, t + At] is then obtained via the inverse Fourier transform 

p n\t) = h I Tr{e £ M A ^ mx e £ (°)Vo}dx- We have used this expression to 

— 7T 

evaluate the fate of the Dicke super-radiance flash ^ (n(t)) at finite time res- 
olution At (not shown). For At > (where V is the single atom emission 
rate), the detector will essentially integrate over the radiation flash such that 
the signatures of collective emission can no longer be resolved in time. In the 
following, we will argue that this situation drastically improves in steady-state 
super-radiance (or more generally, stationary collective transport). 
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3.2 Stationary transport and 'super-transmittance' 

We now consider the case of stationary photon transport between N a = 2 photon 
reservoirs a £ {S,D} (source and drain). First of all, even when the reservoirs 
have non-thermal occupations ng > no, the stationary state of the medium 

P M T z 

(defined by £(0)p — 0) is a thermal state, p oc e ~ : with a temperature 
j3~^ that corresponds to coupling to a single fictitious (N a = 1) reservoir with 
occupation 

T s n s + T D n D _ 1 

n M = =; — -p = -s-jj • 4) 

I s + 1 d e PMil — 1 

Such effective thermalization relations hold quite generally for effective tridi- 
agonal rate equations when the system only supports a single transition fre- 
quency [20]. Specifically, in the high-temperature limit (3 a Q <C 1, this reduces 
to Tm = {TsTs + TdTd)/(J^s + Td), consistent with previous findings for cou- 
pled harmonic oscillators [21]. 

Taking advantage of Eq.(4), the stationary photon current In from the 
source bath through the coherent iV-level medium into the drain bath may 
be evaluated conveniently in the angular momentum basis (p m = (Jm \ p \jm}), 
where one has p m +i = Pm^M I {nM + 1)- From the moment-generating function 
of the drain we obtain (note that trace conservation implies Tr {£(0)^4} = for 
arbitrary operators A) 



N ee lim A (n(t))= (-i9 x )Tr{£( X )e £ WV} = (-i)Tr {£'(0)p} 
= T D (l + n D )Tr{j-pJ+} -T D n D Tr{j+pJ-} 
= (n s - n D ) Fsr ^ <J N , (5) 



with 
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This is one of our central results. First, for N — 1 we obtain 

a \ = r ,,.„ ^tr 11 /! — t , coinciding with e.g. Ref. [22], which in the linear re- 

1 s(l+2ns)+l u(l+2n D ) ' » o L J 

sponsc limit Sn = ns — nry — > with a vacuum drain reservoir (no =0) reduces 
to unity. In this limit, we find the linear 'photon-conductance' lim = 7 c i, 

where 

7d = p , r (7) 

is the classical transfer rate between reservoirs coupled in series via a medium 
with a single excitation. 

Second, for N > 1 and finite ns,nr) > 0, the dimcnsionless quantity <7jv 
describes the degree of collectivity of the steady-state transport process between 
5 and _D. In the large ./V-limit, ojv is asymptotically given by a single parameter 
scaling form, 
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from where the main limits of the current are easily accessible: when the occu- 
pation parameter Um (Eq.(4)) is much smaller than the effective system size N 
(at constant density of two- level systems in the medium), i.e. for um/N < 1, wc 
have a ~ N. In this limit, the current In is an incoherent sum of independent 
contributions from N classical transfer processes with rates 7 c i and thus scales 
linearly with N. 

In contrast, for um/N ;g> 1 the reservoirs provide sufficiently large photon 
flux to drive the medium into the regime of collective, steady-state photon 
'super-transmittance'. There, one obtains to leading order un ps , such 
that the current reads 

n s - n D N 2 

In ~ 7d (9) 

n M 6 

with the quadratic dependence in N indicating the collective nature of the trans- 
mission process. Note that when njj —> (and T$ > 0), this further reduces 
to In « TjjN 2 /2, which gives the theoretical maximum of the current in the 
collective transport regime (the independence of the source tunneling rate Fg is 
a consequence of the required limit of ng — > oo here). Super-transmittance also 
occurs for transport between two thermal reservoirs at sufficiently large temper- 
atures T D = T and T s = T + AT > h£tjk B . For AT -> 0, the classical thermal 
conductance fc t h = liniAT-yo QIn/^-T through the medium is then obtained 
from Eq.(5) and again scales quadratically with N, kth = 7ci (y) N(N + 2)/6, 

whereas at small temperatures T <§C Ml/ks one has k t h = 7ci 

For N > 4, the long-term limit of the cumulant-gcnerating function for the 
current (given by the dominant eigenvalue of C(x) as x ~ > 0) cannot be deter- 
mined analytically anymore. To analyze its global analytic features we therefore 
consider the characteristic polynomial |£(x) — Al| = f(x)i f° r which wc find the 
symmetry /(%) = / (~x - iln ( ^f(jj^ )j- Tnis analytic symmetry transfers 
to the dominant eigenvalue, such that we may deduce for thermal baths the 
steady state fluctuation theorem [23] 

lim Jjl^L = e Ws-p D )» (10) 

t^ooP_ n (t) V ' 

for arbitrary system size N, linking the probabilities of n tunneled particles 
after time t with both reservoir temperatures. 



3.3 Higher photon current cumulants 

We now discuss our second key result concerning the non-trivial scaling behavior 
with N of the higher stationary photon current cumulants Ck = lim ((n fc (f))) 

at k > 1 in regimes where the stationary photon current In = C\ itself is less 
suitable to indicate collective features. 

For large times, particle cumulants ((^"(i))) grow linearly, such that they 
are related to the current cumulants via ((n k (t))) C k t + S k with a constant 
shift term S k ■ Already for moderate system sizes the cumulant-gcnerating func- 
tion cannot be determined directly anymore, calculating the Laplace transform 

of the moment-generating function M(x, z) = Tr |[z • 1 — C(x)] 1 p| however 
only requires the inversion of an (7V+1) x (7V+1) matrix. Therefore, we make use 
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of the relation between cumulants and moments: As the n-th moment is given 
by a polynomial of order n in the first n cumulants, we obtain e.g. for the third 
cumulant ((n 3 {t))) = (n 3 (t)) - 3((n 1 (t))) ((n 2 {t))) - ((n 1 ^))) 3 , which in the 



long-term limit reduces to C3 — > lim ^ 



t— f 00 



(n 3 (t)) - 3 (Cit + 5i) (C 2 t + S 2 ) - {C\t + Si) 



Using limit properties of the Laplace transform, we may therefore directly ex- 
tract the current cumulants C'k once all smaller cumulants and the Laplace 
transform of the fc-th moment is known. 

The most obvious regime is Sn = 0, where there is no net stationary transport 
through the system. In fact, we find that all odd cumulants Cafc+i vanish. Still, 
the two reservoirs (with environmental photon occupations ns = no = tie 
at the systems transition frequency f2) induce stationary fluctuations (thermal 
fluctuations for thermal occupation numbers tie)- Explicitly, for small ue and 
symmetric coupling L5 = To = T for simplicity we find 



C 2 = (1 + r ^ + [N + (N + 2)(N - l)n E + 0{n|}] , 



[(l + nfi)^ 1 

(11) 

where the quadratic scaling corrections oc N 2 underline the collective effect of 
the fluctuations. In the opposite case, where tie is large, the fluctuations behave 
asymptotically as 

C 2 -> ^N(N + 2), 
6 

C 4 -> ^|7V(7V + 2) [iV(iV + 2) + 12] , (12) 

i.e., the super- linear scaling of higher cumulants in the transient Dicke effect 
transfers to the stationary regime. In essence, this already demonstrates that 
even in regimes, where the current does not reveal any collective signatures, 
higher cumulants may do so. 

Wc have also derived analytical expressions for the cumulants in the cases 
of non- vanishing but small and very large bias Sn = ns — Tin, with numerical 
results for finite 5n and njj = displayed in Figure 2. For vanishing drain 
occupation no = and small bias Sn = ns, the cumulant- generating function 
scales as C(x, t) = 7ci(e lx — l)Nnst, which yields a Poissonian process with 
cumulants Ck = JciNns- In contrast, in the large bias limit Sn —> 00, the FCS 
depends quadratically on also for finite no, and we recover the collective 
transport regime: The corresponding cumulant-gcncrating function then reads 

C( X ,t) = ^-N(N + 2)[(e + ^-l)(l + n D ) + (e-^-l)n D ]t, (13) 

which yields the photon current cumulants Ck = TjjN(N + 2)/6 [1 + 2noSk, e ven] 
and describes the independent superposition of two counter-propagating Pois- 
sonian processes (with identical cumulants) corresponding to absorption and 
emission. One recognizes in Figure 2 that indeed higher cumulants yield an ear- 
lier onset of super-linear scaling behavior with the system size N, as compared 
with the stationary photon current In — C\. 
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4 Experimental verification 



Finally, we propose an experimental setup for testing our predictions: In order 
to achieve the small sample limit of super-radiance, it is necessary to constrain 
the size of the optically active medium to dimensions smaller than the photon 
wavelength. We therefore propose a gas of N atoms confined in a cylindrically 
shaped nano-cell [24], where N is controlled, e.g., by evaporation from a particle 
reservoir. We recall that the cross-over between collective (In oc TV 2 ) and 
non-collective (In oc N) photon current regimes is (within the limits of our 
model Hamiltonian) solely controlled by a single parameter, um in Eq.(4). For 
thermal sources, to achieve the collective regime where hm 3> N would require 
large source occupations ns(£l) and thus very large temperatures fceT 3> Ml. 
Therefore, in order to reach this regime we suggest as source an optical cavity 
with loss rate k and pumped by a laser at frequency ojl with pump strength p, 
compare Figure 1 (b). The corresponding stationary cavity photon distribution 

2 

is = ( Q _^yi +K 2 ■ Changing e.g. the detuning of the laser then enables 

one to change ns(£l) and thereby to switch between the different scaling regimes 
discussed above. The nano-cell is placed within the cavity, and in directions 
perpendicular to the cell axis, the atoms are then collectively driven by the 
cavity modes and collectively emit into (vacuum) non-cavity modes which have 
ri£>(f2) w at optical transition frequencies. The photon emissions are counted 
and time-resolved by a detector whose signal is then used to construct the 
distribution P n (t). Sampling a longer signal over different initial times then 
allows one to efficiently construct also higher cumulants Ck [25]. 

5 Summary 

To conclude, we have demonstrated that the transient Dicke super-radiance 
may be transferred to a stationary regime. There, it is reflected as a super- 
linear scaling of transport cumulants with the system size, which we have termed 
super-transmittance. Both in transient and stationary setups, higher cumulants 
yield a much more sensitive signature of collective effects than the first cumulant. 
The effects should be measurable with present techniques. 
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Figure 1: (Color Online) a) Excitation transfers between a medium of N non- 
interacting two-level systems and N a = 2 bosonic reservoirs (source S and drain 
D with occupation ns > no) occur at rates T$/d an d are tracked by counting 
devices, b) Laser-driven optical multi-mode cavity with pump (loss) rate p(k) 
containing ./V atoms that collectively radiate into non-cavity modes. Nanomctcr- 
sized gas cells constrain the atoms in two dimensions and allow for a variation 
of the atom number N in the interaction region. 




Figure 2: (Color Online) Logarithmic plot of the first few current cumulants C\ 
(black), C2 (red), C3 (green), and C4 (blue) in units of T = Ts = ^d versus the 
source occupation ns for different system sizes N = 1 (solid), N = 2 (dashed), 
N = 4 (dash-dotted) and N = 8 (dotted) . The first cumulant shows a relatively 
smooth transition from linear (small differences for small ns) to quadratic (large 
differences for large ns) scaling with the system size N. The higher cumulants 
however show super-linear scaling behavior in a regime where the current does 
not (compare the inset taken at the dashed line ns = 1). Other parameters: 
T s = T D = T, n D = 0. 
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